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Based on N-body simulations, we show that reaUstic galactic disks are subject to bending insta- 
bilities of fire-hose type when the disks are substantially self-gravitating, that is, if they contain 
dark matter distributed in the disk. Depending on the degree of instability, S and U-shaped, as 
well as asymmetric warps are generated. In some cases, the warp may last several galactic rota- 
tions, particularly when the instability is marginal. Since the bending instability is very sensitive 
to the disk flattening, the fractions of dark matter distributed in the disk and in the dark halo are 
constrained. For a Milky Way like galaxy the extended dark halo can not exceed 30-40% of the 
total mass within 35 kpc if the Milky Way warp results from a bending instability. This mode of 
warping provides a unified picture of spiral galaxies, where bars, spiral arms and warps result all 
from disk gravitational instabilities, radial or transverse, which are constantly regenerated by the 
dissipative gas component. 
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Figure 1: Ratio c^/gr as a function of the vertical dispersion at = 15kpc. The values are taken at the 
radius where c^/cr is minimum. The dotted line corresponds to Araki's limit. 



1. Introduction 

Since the first observation of the Milky Way warp at the end of the fifties, warped galaxies 
have represented a challenge for astrophysicists (see for example Binney, 1992, for a review). 
Several hypotheses like normal modes, accretion events, or dynamical friction have been proposed 
to explain their ubiquity. Most of the explanations assume that galactic disks are embedded in a 
dominant spheroidal dark halo. However, there are now several observations (e.g. of NGC 2915) 
which can be best explained with a substantial amount of dark matter in the outer galactic disks 
(Masset & Bureau 2003), in agreement with cold and clumpy gas dark matter models distributed 
like HI disks (Pfenniger & Combes, 1994). In this paper, we summarize the recent result presented 
in Revaz & Pfenniger (2004) where warps are spontaneously generated by bending instabilities 
(see for example Fridman & Polyachenko, 1984, for a formal description of this instability) if a 
substantial fraction of the dark matter resides in the disk of spiral galaxies. 

2. Initial N-body models 

In order to study the possible spontaneous vertical bending of spiral galaxies, we have first 
studied the evolution of equilibrium N-body models starting with different thicknesses, and with- 
out any dark halo. Each model contains a bulge, an exponential disk and a dark matter heavy disk 
parametrized by its thickness h^o. The initial vertical velocity dispersion is obtained by satis- 
fying the equilibrium solution of the Jeans equation, separately for each component /. The radial 
and azimuthal velocities dispersions Or and are obtained from a^, using the ratio between the 
epicycles frequencies K and V, and the rotation frequency Q.: 




(2.1) 



The free factor (3 is chosen of order 1 to keep the Savronov-Toomre radial stability parameter Q 
around 1.5. Fig. [T] shows the Araki ratio o^/or (Araki, 1985) as a function of the vertical dispersion 
a^. Thin disks lie at lower left end while thick disks at the upper right end. 
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Figure 2: Edge-on projections of the models 01, 02, 04, and 06. Time in Myr is indicated at the upper left. 
The box dimensions are 100x25kpc^ 



3. Simulations and results 

As expected from theory, the stability of the models depends strongly on the precise thickness 
of the heavy disk. In Fig. [T], from left to right, models can be divided into four groups. A charac- 
teristic snapshot of each group is displayed in Fig. ^ 

1) Very thin models like model 01 (h^o < 150pc) have a ratio o^/or below 0.2. They are highly 
unstable and generate an asymmetric warp that extends up to z = 4kpc at = 35kpc. 

2) Models 02, 07 and 03 (/z^o = 150 - 250pc) have still a ratio Gz/gr well below the Araki limit. 
The bending instability occurs during the first 2Gyr. A spectacular axisymmetric bowl mode 
(m = 0) grows during about 1 Gyr, before that increases, which stabilizes the disk. 

3) In the models 08, 04, 09 and 05 (/z^o = 250 — 450pc), the ratio o^/gr is almost critical with 
respect to Araki's limit. These models develop S-shaped warped modes (m = 1). In the case of 
model 08, the warp is long-lived and lasts more than 5.5 Gyr. 

4) The thickest model, 06 (/z^o ^ 550pc) has a ratio g^/gr well above 0.3. The disk remains stable. 



4. The influence of an extended dark halo 

The previous models have been run without any outer dark halo, showing that the initial g^/gr 
ratio in the initial N-body runs does determine the warp growth. The influence of an external dark 
halo on Araki's limit is then studied semi-analytically by embedding the previous models in a 
Miyamoto-Nagai potential. The heavy disk to total dark mass ratio / and the halo flattening C^h 
(iso-density axis ratio) characterize the semi-analytical models. For a realistic fixed disk thickness 
hzo = 250 pc. Fig. |3| shows the ratio g^/gr at its minimum value along the radius, indicating the 
stability of the disk. In the absence of the halo (/ = 1), the model is unstable. At the opposite, 
when the heavy disk is replaced by a spherical halo (^/^ = 1) of equal mass (/ = 0), the disk is 
stabilized, as long as the halo flattening is realistic {t^h > 0.5). For very small value of ^z^, the halo 
plays the role of the disk and the model becomes unstable. For a realistic ^, the halo flattening has 
little effect on the stability. In contrast the dark matter fraction in the halo is well constrained: the 
disk remains stable as long as the heavy disk to dark matter mass ratio is smaller than 0.5, in other 
words, bending instabilities appear only if the dark disk is heavier than the dark halo. 
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Figure 3: Stability of a realistic model h^o = 250 pc against bending oscillations. The radial minimum value 
of Cr/Cz is plotted as a function of the dark matter fraction / in the heavy disk and of the halo flattening C^h- 
The solid line corresponds to the Araki critical value of 0.293, the upper dashed line to the values of 0.4, the 
bottom dashed line to 0.2. The shaded region corresponds to unstable models. 



5. Conclusions 

Our numerical N-body simulations have shown that heavy disks may be subject to bending 
instabilities yielding S, U or asymmetric warps. The shape and amplitude of these warps are similar 
to the observed optical warps. This scenario gives also an interesting constraint on the halo mass 
in spiral galaxies. For realistic disk thickness, bending instabilities may appear only if the halo is 
less massive than about 0.3 — 0.4 times the total galactic mass within the warped disk radius. 

As discussed in more detail in Pfenniger & Revaz (2004), the dissipational behavior of the 
gas (not taken into account in the present simulations) reduces the gas velocity perpendicular to 
the galactic plane while, in the plane, higher values of velocities are maintained by the dynamical 
heating due to instabilities like bars and spirals. As a consequence, we expect that a galactic disk 
naturally maintains a regime marginally unstable with respect to bending instabilities. Like spiral 
arms, warps can be repeatedly excited as long as the dissipative component, the gas, exists. 
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